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,-H We estimate large-scale curvature perturbations from isocurvature fluctuations in the waterfall 

(*^ field during hybrid inflation, in addition to the usual inflaton field perturbations. The tachyonic 

^N) instability at the end of infiation leads to an explosive growth of super-Hubble scale perturbations, 

I but they retain the steep blue spectrum characteristic of vacuum fluctuations in a massive field 

^ during inflation. The power spectrum thus peaks around the Hubble-horizon scale at the end of 

^~5 inflation. We extend the usual 5N formalism to include the essential role of these small fluctuations 

f~^ when estimating the large-scale curvature perturbation. The resulting curvature perturbation due 

^SJ to fluctuations in the waterfall field is second-order and the spectrum is expected to be of order 

10~^* on cosmological scales. 

o 

^ I. INTRODUCTION 

O I Cosmological inflation driven by a single scalar field has many attractively simple properties. In particular the 

1^ quantum fluctuations of a slow-rolling scalar field have an almost scale invariant povifer spectrum and these give 

rise to adiabatic curvature perturbations on super-Hubble scales which have an effectively Gaussian distribution for 

canonical scalar fields. This is consistent with current observations of primordial curvature perturbations, but there 

^ is much interest in possible deviations from this simple paradigm which would require either more than one relevant 

'— ' field during inflation or a non-canonical field, or both. Multiple fields can support non-adiabatic perturbations on 

_. 1 super-Hubble scales which can affect the resulting primordial power spectrum and lead to non-Gaussianity. 

►^ Hybrid inflation, driven by the energy density of a false vacuum state which is destabilised when a slow-rolling field 

^f\ reaches a critical value, was originally proposed by Linde [P and subsequently analysed by Linde f? and many others 

l/~j [3 H] • It has proved to be particularly successful for realising inflation in supersymmetric models of particle physics 

^^ [IHZl- In hybrid models there is another field, usually called the waterfall field, which is trapped in a false vacuum 

'^ state until the instability is triggered. Nonetheless the vacuum fluctuations of this field are usually neglected since 

\f^ the field is massive when scales corresponding to large-scale structure of our observable universe leave the horizon. 

^^ Quantum fluctuations of a massive scalar fleld, with mass much larger than the Hubble scale, remain overdamped even 

^D on super-Hubble scales, corresponding to decaying, oscillating mode functions and thus a steep blue power spectrum. 

'^ However it was recently suggested [3 |9] that non-adiabatic large-scale perturbations in the waterfall field could 

^ play an important role, in particular leading to non-Gaussian curvature perturbations ^, since these long wavelength 

modes experience the most rapid growth due to the tachyonic instability in the waterfall field at the end of inflation. 

It is known that a slow transition during hybrid inflation, allowing inflation to continue for some period after the 

jrt tachyonic instability, could lead to large curvature perturbations on scales which leave the horizon around the time 

of the instability iIUl 111] . But it is difficult to model the large-scale primordial curvature perturbations through the 

phase transition where perturbations about the classical background necessarily become large. 

We reconsider this issue in this paper, using a JTV-formalism fl4H17j to evaluate the primordial curvature pertur- 
bation on large scales. We stress the essential role of small, Hubble-scale field perturbations at the end of inflation in 
determining the local integrated expansion N = J H dt in parts of the universe with different values of the waterfall 
field averaged on large, super-Hubble scales. We show that the variance of the waterfall field becomes dominated 
by Hubble-scale perturbations when the tachyonic instability begins, rapidly leading to the end of infiation. The 
duration of inflation is shown to be independent of the large-scale field at first-order, simply due to the symmetry 
of the potential. The curvature perturbation due to long-wavelength modes of the waterfall field are shown to be 
suppressed due to the steep blue spectrum of the waterfall field fluctuations, similar to the case of false vacuum 
inflation supported by thermal corrections |18j . 
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II. WATERFALL FIELD PERTURBATIONS 



We will consider the original hybrid inflation model [H [2l |4] which is described by a slowly rolling inflaton field, </>, 
and the waterfall field, Xi with a potential energy 

y (^, x) = (m' - ^x') ' + Im^^' + h^\' (1) 
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The first term in the potential is a Mexican hat potential for the waterfall field with the false vacuum at x = and 
true vacuum at x^ = 2Af ^/-s/A when </> = 0. The effective mass of the waterfall field in the false vacuum state is 

m2(0)=7(02_02) . (2) 

where we define ^^ = 2\^^^M^ /^. Thus the false vacuum is stabilised for (jp > (p^, while for (jp < (jr^ there is a 
tachyonic instability. 

Note that the simple potential (llj) with a real scalar field, x, has two discrete minima at x = ±2AI/y/X. Thus regions 
which settle into different true vacuum states are separated by domain walls at late times. However vacuum states with 
higher-dimensional vacuum manifolds may have cosmic strings, monopoles or no topologically stable defects. We will 
neglect the formation of cosmic defects while noting their presence could have important cosmological consequences 
in particular hybrid models [H [11] . 

A. Background solution 

In a spatially-fiat Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmology the evolution equations for the back- 
ground fields are 

4> + 3H4> = -(m2+7x')(/. (3) 

X + 3HX - (2a/AM2 - 702 _ Ax2)x (4) 

where the Fridmann equation gives the Hubble rate 

H'^^-^(v{c^,x) + l<P' + lx') (5) 

Initially we assume (/)> (jjc so that the x field is held in the false vacuum and we have the background solution 
X = 0. For simplicity we will work in the vacuum-dominated regime [4] where we can neglect the energy density of 
the (j) field upon the Hubble expansion, H ^constant, 

a ^ H-^ exp [He (t - Q] where H^ = . (6) 

O 

Then the late-time solution for Eq. (Isl) yields 

= 0cexphrff,(i-t,)] (7) 

with 
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'-2-^-4-H-r («) 

The (j) field slow-rolls for m <^ H until the critical point, t = tc, when the tachyonic instability is triggered. In our 
numerical examples we choose r = 0.1. One can verify that this is consistent with vacuum-dominated regime ([6]) so 
long as 7 is not extremely small [1]. 

B. Linear perturbations 

Quantum fiuctuations in the slow-rolling (/> field (coupled to scalar metric perturbations) lead to an almost scale- 
invariant spectrum of adiabatic curvature perturbations on super-Hubble scales during infiation '19j . These correspond 



to local perturbations in the evolution along the classical x = background solution. They give rise to effectively 
Gaussian primordial curvature perturbations usually considered in hybrid inflation models [4]. By contrast quantum 
fluctuations in the waterfall field correspond to isocurvature field perturbations during inflation decoupled from both 
0-field and metric perturbations at linear order |12j . 

Linear ^-field perturbations, with comoving wavenumber k, obey the evolution equation 

5xk + SHSxk +(^~PHl+ ^A Sxk = . (9) 

where the bare tachyonic mass of the waterfall field relative to the inflationary Hubble scale is given by 

P = 2VA— . (10) 

We assume /3 ^ 1 so that the time-scale associated with the tachyonic instabilty is much less that a Hubble time and 
inflation ends soon after the instability begins. If /3 is of order unity or less then there is the possibility of an extended 
period of slow-roll inflation and associated large metric perturbations on scales associated with the transition [TOl E] . 
In our numerical solutions we choose /3 = 100. 

If we substitute the solution for (Eq.Q), and rewrite the evolution equation ^ in term of the canonically 
quantised variable, u — aSx, then we obtain a wave equation with time-dependent mass 

ut+[k^ + ^x\7^)]uk^0 where /.^(r;) ^ ffl^^-^^1^ . (11) 



where primes denote derivatives with respect to conformal time, r] — J dt/a. Note that in de Sitter (l6|) the conformal 
time is given by 

V= F^ = - exp [~H, {t - Q] . (12) 

We have chosen the normalisation for the scale factor in Eq. (l6| such that aH = 1 and 77 = —1 when t — tc, i.e., 
Fourier modes with fc = 1 leave the Hubble horizon when cj) — (j)c and the tachyonic instability is triggered. 
For k — wc find an analytic solution jTT] 



"0 - i-v)'^' 

where r is defined in Eq. (|8| and 
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(13) 



(14) 



In Figure nl we show a typical solution for uo(^)- 

For all modes with a finite k the time-dependent mass-squared in Eq. (Ill is positive, ^J?'{'q) — > /3|7y|^^'^^^^\ but 
negligible at early times where krj — >■ —00. We thus normalise the mode functions at early time to the quantum 
vacuum for a free field, u^ oc e^''^''/v2fc. The time-dependent mass-squared reaches a maximum 

2 (/3 + 2)r //3(l-r)V/'' ^ ^ ( P + iV'^"- ^ ,^, 

and then tends to minus infinity at late times ^^{t}) — > — (/3 + 2)/rf as 77 — >■ 0. At late times we have the asymptotic 
solution 

Uk (X (-r/)-^-i (16) 

corresponding to the tachyonic growing mode solution x = XcGxp [sH{t — tc)] with 



' = \l\+P-l- (17) 
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FIG. 1: The amplitude of the infinite wavelength solution, uo given in Eq. ( |13[ ), plotted as a function of conformal time, ri, for 
r = 0.1 and /? = 100. 



The amplitudes of the mode functions, |u^(77)| for different values of k are shown in Figured 
The power spectrum is defined as 
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T> — \y 



(18) 



Modes with fc ^ 1 leave the Hubble- horizon well before the instability, when \ri\ ^ 1. At these early times the time- 
dependent mass term in Eq. (Ill is positive with /i^(?7) — /3\ri\^'^^^~^' 3> fc^, so that the mode function is suppressed 

12(1- ■ • 



on super-Hubble scales and we find \u \ ex |r/| 
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' . Normalising the initial state to the quantum vacuum 



for a free field then leads to a steep blue spectrum on super-Hubble scales with V^ oc k^, which is even steeper 
than the quantum vacuum for a free field on the smallest scales, V-^ oc k'^. There are no classical perturbations on 
super-Hubble scales before the transition. 

Modes with k ^ 1 remain within the Hubble horizon until after the tachyonic instability is triggered. Even after 
(j) — (pc, gradient terms stabilise short-wavelength modes with fc^ _(_ ^"^(^rj) > and longer- wavelength modes begin to 
grow first. As a result we find that the power spectrum for the waterfall field begins to peak on scales of order the 
Hubble horizon at the transition, A; ~ 1. This is clearly shown in Figure |3j where we fix /3 = 100 and r = 0.1. 



III. APPLICATION OF SN FORMALISM 



In order to calculate the primordial curvature perturbation, (, due to fluctuations in the waterfall fleld we can 
calculate the perturbed expansion, C = 6N = S{hia) |13H17j . from an initial spatially-flat hypersurface to a final 
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FIG. 2: The amplitude of the mode functions, jit||, normalized with respect to their final value, as a function of conformal 
time, rj, for r = 0.1 and /3 = 100. fc = 1 correspo nds to modes leaving the Hubble-horizon when (f> — (f>c. The black dashed line 
indicates the expected late time behaviour, Eq. ( 16l, as (/> — > 0. 



uniform-density hypersurface 
standard first-order result 



For adiabatic perturbations in the (/)-field during the slow-roll phase this yields the 



where HS(j)/(f> can be evaluated on spatially-flat hypersurfaces any time after Hubble exit for each Fourier mode. 
Perturbations in the waterfall field are non-adiabatic field perturbations and hence we must follow their effect on the 
expansion through the tachyonic transition from false to true vacuum at the end of inflation. At first order we can 
simply add these independent contributions to the primordial curvature perturbation from both fields 



c = a + Cx 



(20) 



The usual "separate universe" approach \W is to calculate the classical expansion for different initial values of the 
field assuming a locally homogeneous and isotropic (i.e., FLRW) cosmology, assuming the long- wavelength behaviour 
is independent of much shorter wavelength modes. However in the hybrid inflation model the classical solution 
assuming homogeneous fields is liable to give an incorrect estimate of the duration of the transition. Indeed the 
classical background trajectory assumed in section iHlwas x = for which, classically, inflation never ends and 
N = ln(a) -> oo as t -> oo and cj) -^ in Eqs. ^ and^TJ). 

The dashed line in Figure H shows the classical expansion, Nf{x*) — N{x* — >■ Xf) from an initial value of the 
waterfall field close to the transition, x*j to a given final value, x/, for our chosen values of r and /3 using the 
linearised equations of motion. In particular this shows the singular behaviour of the classical solution for Nf{x*) 
about X* = 0. 




FIG. 3: The power spectrum of Uk as a function of the comoving wavenumber, k, relative to the Hubble-horizon at the transition. 
We show the power spectrum at four different times to show its time evolution. On large scales, k <^1 the power is suppressed 
before the transition, rj < rj,, but then grows rapidly due to the tachyonic instability, for rj > rj^, where rj, ~ — 1. 

Small scale quantum fluctuations play an essential part in the transition as they are amplified by the tachyonic 
transition. The variance of the x field averaged on some scale, L — 27r//c/, > (aH)^^ is given by integrating over all 
longer wavelengths 
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(21) 



On super-Hubble scales we have a steep blue spectrum Vu{k) ex k^ and hence {Sx'^)l = 'Pu{kL)/ia^ ■ We see clearly 
from figure [3] that the variance peaks on scales of order the Hubble scale at the transition. Thus we need to include 
these Hubble-scale modes in our estimate of the background expansion, (N) l, and the perturbation, C = S{N)l, due 
to longer wavelength modes. 

We can always split a Gaussian random field, such as the initial x-field fluctuations at the transition, into a 
long- wavelength part and a short-wavelength part 29 : 



X*ix) ^XLix) + xs{x) = 



f^-in 



d'kxke-'''^ + 



(fkxke~ 



"^sn 



(22) 



In this case we will identify ksput with the peak of the spectrum, around Hubble scale at the transition, /c*. In 
a Gaussian field the long and short wavelengths modes are uncorrelated. Thus we may think of the Hubble-scale 
fluctuations as a statistically homogeneous distribution regularising the divergence of the classical solution for the 
long wavelength field N{xl), as illustrated in Figure |4] 
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FIG. 4: The integrated expansion from an initial time rj, — —1 to a final density corresponding to Xf = W^Hc as a function 
of the average value of the waterfall field, xl, averaged on some super-Hubble scale, L ^ |r;»|. The upper dashed line shows 
the classical solution to the linearised equation (|9| in the long- wavelength limit, k — 0, i.e., assuming a homogeneous field 
X* = XL- The lower dashed line shows the classical solution to the non-linear FRW equations ( 3|5 l. The solid lines shows the 
average expansion in each case integrated over a Gaussian distribution for the Hubble-scale field, x*i with average value xl 
and variance ctJ ~ 0.02H^. Both lines are obtained using r — O.l and /3 — 100. 



A precise calculation of the expansion through the transition would require a full non-linear numerical simulation of 
the inhomogeneous quantum fields on a lattice. We will instead make two important simplifications which will enable 
us to make a simple numerical estimate and which we expect to capture the essential physics. Firstly we will treat 
the fiuctuations as classical on Hubble scales shortly after the transition. This becomes valid due to the tachyonic 
growth of the mode function, but is only marginal at the transition. And secondly we will use the separate universe 
assumption, neglecting spatial gradients on these scales, which again is only marginal at the transition. Modes leaving 
the Hubble horizon at the transition will continue to grow relative to the Hubble scale as inflation continues, but 
eventually will return within the Hubble scale as inflation ends. The separate universe assumption on these scales is 
only marginally valid and for a limited period, but this is precisely the period we wish to study. 

Using these assumptions we can then estimate the expansion, Nf{x*), in a Hubble-scale patch from an initial value, 
X* , around the time of the transition. However in a much larger super-Hubble scale region with an average value for 
the waterfall field, xlj we sample many Hubble-scale patches and integrating over a large volume can be replaced 
by an integral over a Gaussian distribution for the initial values, x*, given a local background value, xl- This is an 
extension of the usual separate universe picture on large scales which incorporates small scale variance, cr^ — {x's)*^ 
uncorrelated with the long-wavelength field. Thus we obtain 



{Nf)i 



dx*Nf{x*)P{x*\^*-,XL)^ 



(23) 



where P{x*\<^*iXl) is the Gaussian probabihty distribution for the local value of the waterfall field in a Hubble-size 
region at time t* soon after (j) = (j)c, 
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(24) 



given the average value on the larger scale, xl, and the variance of the field on Hubble scales, a^ ~ Pu{k*) / a^ , due 
to smaller scale fluctuations, xs{x)- 

The expression (23) for (Nf)]^ gives us the expansion as a smooth function of the large-scale field xl, and thus the 
curvature perturbation can be expanded about the background solution xl = as a function of the perturbation 6xl 
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where, given the Gaussian probability distribution Eq. (24 1, we have 

d{Nf)L 



dXL 

<P{Nf), 



dx 
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^/2^. 



dx* Nfix*) 4 exp 
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(25) 

(26) 
(27) 



From the symmetry of the system (fTl) under x ~^ ~X '^e have Nf{x*) = ^f{~X*) and hence we immediately see 
from Eq. ([26| that 



d{Nf)i 



dXL 



= 0. 



(28) 



Xt=0 



Thus there is no linear contribution to the primordial curvature perturbation from large-scale fluctuations in the 
waterfall field. The leading order contribution to the primordial curvature perturbation from fluctuations in the 
waterfall field on large scales will be second order, and hence non-Gaussian [51 [H]. 



Linearised solution 



Firstly we obtain numerical solutions for the classical evolution, Nf{x*), from an initial value of the waterfall field, 
X* when cj) — (j)c, to a final value Xf ^ X* using the linearised equation of motion (|9|) for Sxk — X*, where we work 
in the long- wavelength limit and set fc = 0. As expected, the classical evolution Nf{xL), shown by the dashed line in 
figure |4J is singular in the limit x/, — >■ 0. By contrast the solid line shows how the classical singularity at xz, = is 
regularised by the quantum dispersion of the local value, x*j once we include smaller Hubble-scale modes, yielding a 
finite value for {Nf)^ as xl -^ 0. 

We can obtain a rough analytic estimate for Nf{xL) by noting that at sufficiently late times after the transition, 
— ?► and we expect x to have the late-time behaviour, x ^ {~v)^'^ given in Eq. (16 1. Therefore, we expect Nf and 
Xs Ix* to be approximately given by 



A^f = - In ^ + const, 
s X* 



(29) 



Note that this assumes, <C 0^ by which point the linear approximation for x is expected to have broken down. 



Nonetheless, substituting Eq. (29 1 into Eq. (231 allows us to give an analytic estimate 



(A/)i = A^, + ^ln 
2s 



lot 



ln2 + 7EM/2 1 Sx 



x/2s 



2v^S Cr2 



O 



5x 



(30) 



where 7em — 0.57721 is the Euler-Mascheroni constant. 

For our chosen parameters, r — 0.1 and /3 = 100, we have al 
analytic approximation ( 30 1 then yields 



Vulal ~ (Hc/27r)2, and we set Xf = lO^ifc- The 
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const. - 1.62^ 
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(31) 



Performing the integration numerically in the linearised approximation one obtains 



{Nf)L c, const. - 4.62^ + ^ (^ ) • (32) 



So the analytic approximation ( 29 1 , gives a result roughly a factor of 3 away from the numerical result in the linear 
analysis. 

B. Non-linear FLRW solution 

As we have neglected gradient terms by setting A: = in the linear equation for x(-^)j we can in fact solve the 
full non-linear equations in the homogeneous limit, treating the local Hubble-scale patches as separate universes |16| 
obeying the FLRW equations ^^ in order to determine the coupled evolution of x and (f>, and hence N, starting 
from X = X* when = 0c- In this case we must also specify the energy scale of inflation. In Figure |4] we show 
the classical solution Nj{xl) as well as the averaged (Nf)]^ including the small scale dispersion, where we have set 
M ~ 10^^ GeV and hence He — 10^^ GeV. The results are qualitatively similar to those obtained from the linear 
approximation. Non-linearities, in particular the x-dependent effective mass for the (j) field, leads to a slightly more 
rapid transition, especially for larger initial values of xl, and we find 

{Nf)L c, const. - 6.59^ + O (^-^) . (33) 



Note that in our non-linear, separate universe solutions we have calculated Nf up to a final fixed density, and have 
verified that in practice this coincides with a final value for the waterfall field Xf — ^O^Hc- For our chosen parameter 
values the universe is still inflating at this final time, with slow-roll parameter e = —H/H^ ~ 0.001. Thus our initial 
comoving Hubble-scale at ?/* is still larger than the final Hubble scale. But the slow-roll parameter is growing quickly 
and inflation ends soon afterwards (A^^end — Nj + 0.4). 

Independently of the details of the calculation we see that the curvature perturbation due to large-scale perturbations 
in the waterfall field is second-order, and suppressed relative to the Hubble-scale variance of the field at the transition. 

Vc(k)^H-*Vs^2{k). (34) 



where 



n.;w=g/^'^' ^"f'^'g"'^-r" - w 



27r / fc'3 |k-k'|'3 



Given the steep blue spectrum of the super-Hubble perturbations in the waterfall field, Vsxhi^) — 'Pu{k)/a^ oc k^ 
(which implies a white spectrum in the standard terminology since Ps^lW — 27r^7'^^^(fc)/fc3 ex const.), we thus 
conclude that the spectrum of the resulting primordial curvature perturbation, <^^, on super-Hubble scales is also blue 
for fc <C A:* : 

T'c, (fc) ^ H-^Vs^, {K)V5^, (fc) ^ {k/Kf . (36) 

Assuming cosmological scales leave the Hubble-horizon around 40 e-folds before the end of infiation we have 

IV. CONCLUSIONS 

We have estimated the primordial curvature perturbation produced by fluctuations in the waterfall field during 
hybrid inflation. We have calculated linear perturbations about the classical background trajectory, x = 0, during 
slow-roll inflation and then studied how these affect the primordial curvature perturbation when the waterfall field is 
released and infiation comes to an end. To do this we have used a extension of the usual (5A^-formalism that identifies 
the primordial curvature perturbation with the perturbation in the local expansion on a uniform-density hypersurface, 
C^SN. 

The challenge for the standard SN approach is that the homogeneous classical background solution, with smooth 
fields on scales much larger than the Hubble-horizon, fails to provide a good description of the dynamics. A classical 



10 

solution which starts precisely in the false vacuum state stays there and inflation never ends. Thus even a tiny 
perturbation away from this background solution results in a huge apparent change in the local expansion, N. In fact 
the false vacuum is destabilised mainly due to fluctuations of the waterfall field on scales close to the Hubble-horizon 
at the end of inflation. 

The coupling between the slow-rolling inflaton and the waterfall field that leads to the tachyonic instability also 
makes the waterfall field massive before the transition. Thus quantum fluctuations of the waterfall field are not 
amplified during inflation and remain in a quantum state even on super-Hubble scales until the critical value of the 
inflaton is reached. Once the tachyonic instability is triggered, there is an explosive growth of long-wavelength modes, 
but they retain the steep blue spectrum on super-Hubble scales. On the other hand much smaller scales, well inside the 
Hubble horizon at the transition are stabilised by spatial gradients and remain in their vacuum state. The resulting 
power spectrum for the waterfall field thus peaks on scales around the Hubble scale at the transition, as shown in 
Figure [31 and these modes play an essential role in the dynamics that must be included when calculating the effect of 
large scale fluctuations in the waterfall field. 

Thus we identify the primordial curvature perturbation due to perturbations in the waterfall field, Sxl on some 
large scale L, with the perturbation in the average expansion, (Nf) , including fluctuations in the waterfall field on 
Hubble scales at the transition. We adopt a Gaussian distribution for the field whose average value in a region of size 
L is Sxl, but whose variance is given by crj = Vu{ki,) / a? . We find two main results: 

• from the symmetry of the potential (IT]) under x '^ ~X '^s see that the primordial curvature perturbation is 
independent of Sxl at first order, and the curvature perturbation is second-order in Sxl- 

• the spectrum of the primordial curvature perturbation due to fiuctuations in the waterfall field on large scales 
is suppressed by a factor of order {k^,L)^ which for cosmological scales is likely to be of order 10~^^. 

While we have considered the specific hybrid potential ([I]) and presented numerical solutions for specific parameter 
choices we believe the general conclusions will hold for all hybrid models in which the waterfall field is massive during 
slow-roll inflation and for which the end of inflation occurs due to a rapid tachyonic instability. 

It is not surprising that the effect of small scale fluctuations must be taken into account when estimating the 
primordial curvature perturbation when the large-scale field is very close to zero. This is the case, for example, in the 
curvaton model where the curvature perturbation, (^ ~ Sa/a, could become very large when the background field, cr, 
is close to zero in some regions of the universe |22j . The apparent singular behaviour of the curvature perturbation 
in this case is regularised by smaller scale fiuctuations in the curvaton field [53]. 

In our analysis we have considered only the early stages of the tachyonic instability. Eventually inflation ends and 
the coupled flelds oscillate about the true vacuum (/> = and x^ = 2M^/vA. Numerical simulations are required 
to study resonant particle production and other non-perturbative effects in this regime which goes beyond the scope 
of this paper. Resonance is often most efficient in long-wavelength modes, but fields which are massive throughout 
slow-roll infiation necessarily give rise to a steep blue spectrum of perturbations and hence the large-scale power is 
suppressed with respect to smaller scales [24VI26J . 



Note added: While this work was in progress, two papers appeared |271 128] which study the generation of large 
scale curvature perturbations at the end of hybrid inflation from complementary viewpoints and whose conclusions 
appear to be consistent with our results. 
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